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ON AN EXTENSION OF A STATIC BALL SOLUTION 
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A new exact static interior solution of the Einstein equations is obtained for a gravitating ball filled with a Pascal 
perfect fluid . The solution is an extension of the well-known interior solution with a parabolic distribution of mass 
density and describes more compact astrophysical objects such as neutron and hyperon stars. It is shown that a 
behaviour of the new mass density distribution near the stellar centre can be described as a cusp catastrophe. 



1. Introduction 

The problem of finding exact solutions of Einstein's 
equations is connected, above all, with their nonlinear- 
ity. Therefore every new exact solution always calls 
certain interest, both from the point of view of a phys- 
ical interpretation and from the point of view of the 
method of its derivation. 

Schwarzschild's known interior solution (see [1]) is an 
exact interior solution of the Einstein's equations with 
a homogeneous distribution of a Pascal neutral perfect 
fluid. This solution describes the interior field of a static 
star. The main shortcoming of this solution is that the 
velocity of sound exceeds the velocity of light. The mass 
density is constant in the whole volume of the astrophys- 
ical object and has a finite jump on the stellar surface. 

Another exact interior solution [2], [3] is free of the 
above shortcoming. The solution describes a fiuid ball 
with the parabolic law of mass density decrease from 
the centre to the surface. However, both solutions near 
the stellar centre behave almost equally with rrespect 
to Petrov's algebraic type. The Schwarzschild interior 
solution has the algebraic type 0. The second solution, 
in a small neighbourhood of the ball centre, has approx- 
imately the same algebraic type 0. 

In the present article, the problem of generalization 
of the solution with a parabolic distribution of the mass 
density for a neutral Pascal fluid is considered. 



2. Metric and gravitational equa- 
tions 

The Einstein equations with the static metric 

ds^ ^ G^{x)dT^ +2L{x)dTdx~ x^{de^ + sine^d(p^) (1) 

and the energy-momentum tensor of a perfect Pascal 
fiuid 



Tap = (fl+p) UaUp -pgap 



(2) 
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are considered. Thr function /i is the mass density and 
p is the pressure; u" = dx'^ / ds is the fiuid 4D velocity; 
a, (i = 0,1,2,3. 

These equations can be reduced to the nonlinear 
equation of spatial oscillations 



where C, is the new variable 



(3) 



dC = dy/^e, (4) 

with e — G"^ jl?; t — t/R and x — r/R are dimen- 
sionless variables; < 2^ < 1; t,r are the time and 



radial coordinates, y 



< y < 1, and R is the 



stellar radius. The vacuum velocity of light is chosen to 
be equal to unity. 

We can write 57^ in a more convenient form: 



= -d{'^/v)dy 



(5) 



where 



$=!-£ = 



j Ky)Vydy = f / iJ-ix)x'^dx (6) 



plays the role of the Newton gravitational potential in- 
side the ball; x = xi?^. 

When Jl^ equals zero, we have the Schwarzschild 
interior solution. For fi^ = fig = const, Eq.(3) is 



nlG = o, 



(7) 



and we have the known exact static solution [2], [3] 

G ^ y/g^ = Go cos{nQ C{x) + a) (8) 

with flo = X a-nd the parabolic distribution of the 

mass density 

fj. = fiQ{l- a x'^) = 1.10 {I- ay), (9) 

where a = (/io - /i-R)/Aio < 1; Mo = = 0); ^fi = 
fi{x = 1). 

This solution describes compact astrophysical ob- 
jects such as neutron stars. 
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3. Darboux's method and an ex- 
tension of the solution 

According to Darboux's method [4] , an extension of this 
solution is connected with the choice 



with fin = const and 



w{0 = 2{lnYyi^, 



(10) 



(11) 



where F is a particular solution of the equation of a 
linear spatial oscillator (7) with Oq = 6^ = const. Here 
we select the function Y = cos{bQ). 
So we have the extension of Eq.(7) 



G'U + [^0 - 262sec(6C)] G = 0. 



(12) 



The exact general solution of this Darboux equation is 
[5] 

G(C)=A-[r!o -005(00 C + ao)+ 

+b-tg{b(:)-sin{noC + ao)\. (13) 
The new function ^{y) is written as a quadrature: 



Hy) = {xi^omy-^W + b''y J sec^{bx{y))dy. (14) 

With the aid of this function we easily find the new mass 
density distribution 

XAt(y) = XMo - 5ng y + j sec^{bC{y))dy+ 



46\sec2(6C(y)), 

where /i(0) = /Uq. 

Now the new variable ( is the quadrature 



(15) 



C(2/) = 



dy 



-2b'y 



dy ^ 
Vl-%) 

1 - (XMo/3)t/ + Ogy" 
dy 



1/2 



On the surface of the ball we have 
l>(2/ = 1) = r/o + 26^ 



(16) 



(17) 



where % = ^m/R is the compactness for the model 
with a parabolic distribution of mass density; m is the 
integral mass of ball. 

The expression for the new mass density near the 
stellar centre can be rewritten as 



Mo 



125ag - 7846- ^2 



362 



1 



125a^ - 78462 36^ 



(^0 



-y + i^y , 



(18) 



where (tq — x ■ ^J'0■ 

This equation can be reduced to the equation which 
describes the cusp catastrophe 

+ (/c - f)y + Bsi = 0, 

where ei = 1 — fL{y)/no -C is a small deviation from 
the mass density at the stellar centre, fc = 20/3, B = 

(3/2)62//ao. 

The behaviour of the function jl{x) describes the 
star's stability near of the centre under the control pa- 
rameters /xo,6, X- There exists a range of the control 
parameters in which the system has a catastrophic be- 
haviour. 

4. Summary 

Application of Darboux's method to a known exact 
static interior solution of the gravitational equations 
with a parabolic mass density distribution makes it 
possible to obtain a new interior solution with a new 
mass density distribution law. The extension of solution 
can be used for simulation of more compact astrophysi- 
cal objects: neutron stars and, probably, hyperon stars. 
Another feature of the new solution is that near the stel- 
lar centre the mass density behaviour is described by a 
cusp catastrophe, and, for some values of the parame- 
ters included in the cixprcssion for the mass density, the 
stability of the star breaks down. 
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